Insulator-metal transition in the rutile-based VO2 
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The metal-insulator transition (MIT) in paramagnetic V0 2 is studied within LDA+DMFT(IPT), 
which merges the local density approximation (LDA) with dynamical mean field theory (DMFT). 
With a fixed value of the Coulomb U = 5.0 eV, we show how the MIT is understood in a new picture: 
spectral weight transfer accompanying the increase in the displacement of V ion (_L c) within the 
strong correlation scenario. Within this new scenario, we find good quantitative agreement with (i) 
switch of the orbital occupation of (xy,yz + zx,yz — zx), (ii) thermodynamics, and (Hi) the one- 
electron spectral function in the metallic phase of VO2. We also compare our results for the total 
spectral density with other approaches which use QMC to solve the impurity problem of DMFT. 



PACS numbers: 71.28+d,71.30+h,72.10-d 
I. INTRODUCTION 

Electron correlation-driven metal-insulator transitions 
(MIT) in the 3d transition metal (TM) oxides remain 
a problem of enduring interest in condensed matter 
physics. 1 Vanadium oxides have proved to be fascinat- 
ing candidates in this context. 2 Recent X-ray absorption 
(XAS) studies 3 on V2O3 have led to a drastic revision 
of the traditional picture 4 of the MIT in this system, 
forcing a description in terms of a 5 = 1, multiband 
model. 5,6 Similar behavior is also observed in pressur- 
ized Ca2Ru04, 7 suggesting that such a scenario may be 
a manifestation of the general importance of orbital cor- 
relations in TM oxides. 

Stoichiometric vanadium dioxide (VO2) is another ex- 
ample of a system showing a MIT with reduction of tem- 
perature. The high-T metallic (M) phase has rutile (R) 
structure, which changes to monoclinic at the first-order 
MIT occuring around Tmi = 340 K. This, along with 
the observation of spin dimerization along the c-axis, 8 
gave rise to theories linking the MIT to strong electron- 
lattice coupling (ELC) , 9 where the important role of the 
antiferroelectric displacement of the V ions (_L c) was 
pointed out. It was suggested that c-axis dimerization of 
the V-V pairs would open up a gap in the low-T insulat- 
ing (/) state in the framework of a bonding-antibonding 
splitting. The M state, accompanied by disappearance of 
this dimerization, would result from an upward shift of 
the d xy band (d a in Ref. 10). In a classic work, a zero-gap 
insulator was obtained within LDA; 11 this study treated 
structural aspects in full detail. However, the inability 
to reproduce the observed gap value suggested that in- 
clusion of on-site correlations was necessary. In contrast, 
in a pioneering work, the importance of electronic corre- 
lations, along with that of the antiferroelectric displace- 
ment of V ions from their c-axis positions in driving the 
MIT was emphasized. 10 From a detailed perusal of exper- 
imental results known at that time, it was argued that o 
axis spin dimerization is a consequence of the MIT driven 
by electronic correlations in concert with the antiferro- 
electric displacements. Subsequently, a model calcula- 



tion 12 explicitly realized this suggestion, but with unre- 
alistically small model parameters. On the experimental 
side, the observation of anomalous (correlation driven) 
features 13-15 and results from magnetic measurements 16 
finding local moments in Vi_ x Cr x 02, required a correla- 
tion based scenario. Simultaneously, observation of only 
a modest mass enhancement in the R phase seemed to 
suggest that while correlations were undoubtedly present, 
they did not play a crucial role in the MIT. In the light 
of this conflict, the origin of the PI state as well as the 
mechanism of the MIT is an open and controversial issue. 

Recent (unpublished) XAS results 17 show that the 
abrupt MIT in VO2 is accompanied by a hysteresis in 
occupation of different t2 g orbitals, directly confirming 
the importance of orbital correlations in a multi-band 
situation. 5 Recent photoemission measurements show a 
large spectral weight transfer from high- to low ener- 
gies, on a scale of 0(2 eV), across the MIT. Further, 
this character of the MIT, which is driven by a sudden 
increase in carrier density rather than their mobility, is 
shown very clearly by time-delay measurements. 18 More 
recently, Lim et al. 19 found the MIT to be strongly first 
order by I-V measurements on thin films. Very interest- 
ingly, micro- Raman scattering experiment clearly shows 
that this electric field driven Mott transition is not ac- 
companied by any structural transition, in contrast to 
the usual T driven one, strongly suggesting the Mott- 
Hubbard scenario. These observations put particularly 
strong constraints on an acceptable theory of the MIT 
for VO2: (i) it should be able to describe this switch in 
orbital occupation at the MIT, (ii) the metallic phase 
should be moderately correlated (see, however, 15 ), and, 
(Hi) the insulator should be of the Mott-Hubbard type, 
accompanied by c-axis dimerization of V 4+ -V 4+ pairs, 
and the MIT itself should be first-order, consistent with 
Ref. 18. 

In this article, we discuss a concrete theoretical sce- 
nario, extending 10,12 earlier works in the light of the 
above requirements. Using the recently developed ab- 
initio LDA+DMFT scheme, we show how a consistent 
treatment of strong, dynamical multi-orbital correlations 
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along with the antiferroelectric displacements of V in the 
real LDA bandstructure leads to a unified understand- 
ing of the MIT in VO2. We also compare our results 
with other, recent LDA+DMFT work which uses quan- 
tum Monte Carlo (QMC) method to solve the impurity 
problem of DMFT. Finally, we show how very good quan- 
titative agreement with PES, XAS, thermodynamic and 
magnetic measurements is obtained. Our results support 
the Mott-Hubbard picture of the MIT in V0 2 , where the 
Peierls instability arises subsequent to the M-I transition. 



II. MODEL AND SOLUTION 

We start with the LDA bandstructure of VO2 in the 
monoclinic crystal structure corresponding to the insulat- 
ing phase. 20 The LDA DOS, computed using the known 
value of the antiferro-electric displacement of V 4+ ions 
(e, giving the splitting A Qj a between the xy, yz ± zx or- 
bitals at the local level) (Fig. 1) shows interesting fea- 
tures: the lowest lying xy (in the M notation) band is 
the most heavily populated, while the yz ± zx bands are 
less populated. The bonding-antibonding splitting 9 is 
clearly visible in the results, which also show that the 
total bandwidth is about W ~ 2.3 eV, contrary to much 
smaller previous 9 ' 10,12 estimates based on model calcula- 
tions. Except for their inability to reproduce the Mott 
gap, these LDA results agree with Goodenough's argu- 
ments. As emphasized by Mott, 10 it is inconceivable that 
e alone could open up a charge gap (Eq — 0.6 eV), this 
requiring proper incorporation of correlation effects. 




FIG. 1. LDA partial density of states for the d yz ± zx (dot- 
ted) and d xy (solid) orbitals, obtained from Ref. 20. 

Thus, the LDA (bandstructure) part of the Hamilto- 
nian is 



H =J2 e a/3( k ) c L<T C k/3<r + ^e M 
ka/3(7 ia< 



(i) 



where a, [3 = xy,yz + zx,yz — zx. To avoid double 
counting of interactions already treated on the average 
by LDA, we follow 21 to write 



h lda = J2 e M k ) c La c W° + ' 
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where e° ac — - 1 7 ' • 1 

U, Jh as defined below, and the full Hamiltonian is, 
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Constrained LDA calculations give U ~ 4.5 eV, Jh = 
1.0 eV and U' = U a0 = [U - 2J H ) = 2.5 eV. Given 
the uncertainty in the precise estimation of U , we choose 
U = 5.0 eV in the calculations. Considering only the 
tig manifold, the PM/para-orbital phase, and working 
in the basis which diagonalises the LDA density matrix, 
we have G al3aa > (w) = 8 a p8 aa >G aa {<jS) and Y, al3aa >(uo) = 

In the t2 g sub-basis, a DMFT solution involves (?) re- 
placing the lattice model by a multi-orbital, asymmetric 
Anderson impurity model, and (ii) a selfconsistency con- 
dition requiring the impurity propagator to be equal to 
the local (fc-averaged) Green function of the lattice, given 
by 
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(4) 



Using the locality of S Q/ g in d = 00, we get G a (u>) = 
G° (u> — T, a (iv)). Further, since U a p, Jh scatter electrons 
between the xy, yz + zx, yz — zx bands, only the total 
number, nt 2g = J2 a nt 2 g ,a is conserved in a way consis- 
tent with Luttinger's theorem. 

To solve the impurity problem of DMFT, we use the it- 
erated perturbation theory (IPT), generalized to the case 
of t2 g orbitals for arbitrary filling. 22,23 The local propa- 
gators are given by 
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(5) 



Local dynamical self-energies E Q (a>) are computed 
within the extended multi-orbital iterated-perturbation 
theory (MO-IPT) which generalizes the IPT for the one 
orbital case to the multi-orbital situation, for any tem- 
perature 24 and band-filling. 23 Explicitly, 



where, for example, 
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III. NUMERICAL RESULTS 



and G°H = 



The multi-orbital dynami- 



cal bath A a (oj) function is determined using the DMFT 
self-consistent condition which requires that the local in- 
purity Green function to be equal to the local Green func- 
tion of the lattice (Eq. (5)). Finally, in Eq. (6), 



and, 



n a (l - 2n a ) + D ai [n] 
<(!-<) 
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Here, n a and n° are particle numbers determined from 
G a and G° respectively, and D ai [n] ) is calcu- 

lated using 22 



(n a n^) = (n a )(n 7 ) 
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U a jTT 



f(ijj)Im[T, a (uj)G a (uj)]duj . 

(9) 



The last identity in Eq. (9) follows directly from the equa- 
tions of motion for G a (tu). These equations satisfy the 
Luttinger sum rule, and reproduce the first few moments 
of the spectral functions 25 , guaranteeing the correct low- 
and high energy behavior of the propagators. As shown 
by Kakehashi et al. 25 the IPT for arbitrary filling can be 
derived in a controlled way by truncating the memory 
function matrix in a Mori-Zwanzig projection formalism. 
The above equations are solved selfconsistently respect- 
ing these constraints. 

Below, we use the following strategy to study the Mott 
M-I transition. Since we do not expect much correspon- 
dence between changes in bare LDA quantities with those 
affected in non-trivial ways by dynamical electronic cor- 
relations, we choose the bare LDA DOS for the "insu- 
lating" solution (we do not change this DOS to study 
the MIT) at low T(~ 75 K) in our calculation. To ac- 
cess the MIT as a function of T, and driven by change 
of orbital occupation, we monitor the fully renormalized 
spectral functions, their renormalized occupations, along 
with changes in the renormalized antiferroelectric dis- 
placement (A Q( g). To be precise, we start with a trial 
value of e. The multi-orbital correlations have two im- 
portant effects in this situation. Given a trial value of 
e, we expect U a p to renormalize its bare value (the ef- 
fective "field" in the ti g sector) via multi-orbital Hartree 
shifts, leading (t) to changes in the orbital occupations, 
and more importantly, (it) to large changes in spectral 
weight transfer (driven by dynamical nature of strong lo- 
cal correlations) upon small changes in A Q| g, leading to 
stabilization of the second 10 (metallic) DMFT solution. 
In other words, this large SWT drives the abrupt MIT 
around a critical value of A Q( 3 = A c , as is indeed con- 
firmed below. 
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FIG. 2. t 2g partial density of states using LDA+DMFT for 
U = 5.0 eV in the PI (solid) and PM (dashed) phases of V0 2 . 
The upper panel corresponds to d yz ± zx orbitals, and the lower 
panel, to the d xy orbital. 

In Fig. 2 (solid lines), we show our result for the par- 
tial DOS in the insulating case for T — 0.0075 eV. Notice 
that the charge gap, Eq ~ 0.5 eV, consistent with the 
optical data. 26 In accordance with the orbital splitting as- 
signment from LDA+U, the lower-lying xy band is more 
localized in the solid. The orbital occupations in the MI 
state are computed to be (n\\ , — 0.36, 0.32, 0.32), quite 
different from the XAS estimations. From the computed 
spectral functions, the renormalized value of the distor- 
tion (caused by A Q ^) is found to be A^ ~ 0.17 eV. 
Decreasing the trial values of nil (leading to decrease 
in A£g) leads to a large transfer of spectral weight via 
DMFT in the multi-orbital situation, and, at a critical nil , 
to an abrupt (first-order) transition to a metallic phase 
(Fig. 2, dashed lines, computed for T = 0.03 eV). Inter- 
estingly, the d\\ DOS still represents insulating behavior, 
while the d^ DOS exhibits characteristic metallic behav- 
ior, providing an explicit realization of the two-fluid sce- 
nario 2 proposed earlier to describe the general problem 
of MIT. In accord with the above argument, the renor- 
malized A^a ~ —0.29 eV, corresponding to a reduction 
in A a/ 3, in full agreement with Mott. 10 The difference 
6 = A£g — A^g = 0.46, in good agreement with the 
value of 0.5 — 0.6 from cluster calculations. 10 Moreover, 
this agrees with the orbital splitting assignment in the M 
phase, 9,10 ' 20 where the d v band(s) lies lower than the dy 
band. Finally, the orbital occupations in the M phase are 
found to be (niprv) = (0.34,0.33,0.33) respectively, in 
very good semiquantitative agreement with recent XAS 
results. 

Based on these results, the abrupt MIT in VO2 is un- 
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derstood as follows. A decrease in the LDA value of nil 
leads to an effective reduction in A 1 (caused by Hartree 
shifts due to U a p) leading to an effective increase in n v , 
and, more importantly, to large changes in dynamical 
SWT from high- to low energy, driving the first order 
MIT. The converged values of % (obviously different 
from the trial values) vary in a way qualitatively con- 
sistent with indications from recent XAS measurements. 
Thus, in our picture, the MIT is driven by (i) abrupt 
changes in n\\,n v and A Q ^ due to large change in dy- 
namical SWT, and, (it) an increase in the carrier con- 
centration, rather than their mobility, consistent with 
time-delay measurements. 1 
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FIG. 3. (color online) The integrated photoemission line- 
shapes for the PI (solid) and PM (dashed) phases within 
LDA+DMFT for U = 5.0 eV. The LDA result is the 
dot-dashed line. 

Given the relatively wide ti g bands in the LDA (con- 
trast with Ref. 12, where the bandwidths, as well as 
U, U a p, Jh were severely underestimated), correlation in- 
duced mass enhancement is found to be moderate. In- 
deed, from the self-energy (not shown), the average ef- 
fective mass enhancement factor is m* /m ~ 3.0, in full 
accord with that extracted from reflectance studies. 27 
Clearly, this enhanced mass comes solely from the 
band (this is the only band showing "metallic" behav- 
ior). 10 Further, the high-T static spin susceptibility is 
enhanced by the same factor of 3 coming from renormal- 
ization effects of multi-orbital correlations, in contrast to 
the Stoner mechanism discussed before in Ref. 10. These 
results can be interpreted to mean a moderate renormal- 
ization caused by the additional screening provided by 
the wide d v bands in the M phase. 

In Fig. 3, we show the integrated photoemission (PES) 
spectra calculated in both the I and the M phases. Com- 
paring our results with published experimental work, 1 
very good semiquantitative agreement between theory 
and experiment is seen. In this context, we emphasize 



that a certain amount of caution should be exercised 
while comparing theory to experiment because (i) re- 
duced co-ordination at the surface would normally en- 
hance correlation effects, decreasing the DOS at Ep, and 
(ii) the surface of VO2 is known to show insulating fea- 
tures (Ref. 14). Incorporating spin-Peierls effects in the 
PI phase will also change the PES lineshape in the PI 
phase somewhat; this is left for the future. With these 
caveats in mind, the overall agreement is strkingly good, 
and, in particular, the transfer of spectral weight across 
the MIT is faithfully reproduced by theory. 

In particular, in our picture, the small peak in the 
metallic phase is not a quasicoherent Fermi liquid (FL) 
peak, but part of the tail of the d yz ± zx bands (from 
Fig. 2) pulled down below Ep in the PM phase, as de- 
scribed above. More detailed orbital-selective spectro- 
scopic investigations are needed to confirm this picture. 28 
This is also consistent with the absence of a T 2 term in 
the dc resistivity, which is known 15 to follow a linear-in-T 
behavior for 330 K < T < 840 K. We notice that a FL 
peak in PES would contradict the incoherent, linear-in-T 
(non-FL) resistivity. In any case, in correlated systems 
on the verge of a Mott transition, such a FL resonance 
should be visible, if at all, only at very low temperatures, 
and is rapidly washed off with increasing T. 

IV. COMPARISON WITH LDA+DMFT(QMC) 
RESULTS 

We turn now to a comparison of our work with other, 
recent approaches which use QMC to solve the im- 
purity problem of DMFT. 29-31 First, we clearly see 
in Fig. 4 that the QMC calculations resolve a huge 
peak near Ep. Though a miniscule pseudogap can 
be discerned in Ref. 29, a direct comparison between 
both LDA+DMFT(QMC) results with the correspond- 
ing LDA ones clearly shows that the many-body renor- 
malised DOS, p(Ep) at the Fermi surface, coincides with 
the corresponding LDA DOS, even at T — 770 K. (Note 
that the QMC calculations done at T w 500 K 30 are also 
pinned to the same value at EfV). This is clearly prob- 
lematic: let us recall that the Friedel Luttinger sum rule 
implies p{E F ) = plda(Ef) only at T = 32 in d = 00. 
Further, in a DMFT framework, this pinned (at T = 0) 
feature actually ceases to exist above the lattice coher- 
ence scale, T co h, which is quite small near a correlation- 
driven MI transition. In the LDA+DMFT framework, 
such a pinned feature can only be associated with a cor- 
related FL peak. Its persistence up to T = 770 K (QMC) 
is thus in direct conflict with the bad-metallic, linear-in- 
T resistivity observed from much lower T(= 340 K) in 
VO2. Optical measurements further confirm this: 26 no 
coherent Drude peak is observed in the PM phase. 

The pinning observed in LDA+DMFT(QMC) work is 
thus very hard to justify, both on theoretical and ex- 
perimental grounds. Further, explicit consideration of 
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inter-site V-V interactions (dimerisation) within clus- 
ter DMFT would normally cause a low-energy pseu- 
dogap to appear around Ep, making any pinning im- 
possible, even at T = 0. Clearly our results do not 
suffer from this conflict, and are moreover in excellent 
quantitative accord with PES/XAS results, as well as 
with the effective mass enhancement estimated from op- 
tics/thermodynamics, as discussed before. To compare 
various approaches, we show in Fig. 4 our result, together 
with those obtained from LDA+DMFT(QMC) 29 " 31 in 
comparison with recent experimental results. While our 
result quantitatively describes the full spectral function 
from -3.0 < u < 1.2 eV, the LDA+DMFT(QMC) result 
does not capture the detailed lineshape correctly over the 
whole range. 
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FIG. 4. Comparison of theoretical LDA+DMFT(IPT) and 
LDA+DMFT(QMC) result for the total one-electron spectral 
function in the metallic phase of VO2 to the experimental re- 
sults taken from Ref. 33 (for PES) and from Ref. 34 (for XAS). 
The LDA+DMFT(QMC) results are taken from: Ref. 29 (a) 
[T = 770 K], Ref. 31 (6) [T = 770 K], and Ref. 30 (c) 
[T w 500 K]. 

What are the sources of the discrepancy be- 
tween our results and the LDA+DMFT(QMC) ones? 
There are two important differences between our ap- 
proach compared to the LDA-fDMFT(QMC) works: 
The LDA+DMFT(QMC) papers derive the many- 
body density of states (DOS) in the two phases 
(dimeriscd/metallic) by using the corresponding LDA 
DOS separately for each phases, and using DMFT(QMC) 
to study dynamical correlation effects in each phase. 
Thus, the two phases can be studied separately, but it 
is difficult to see how the first-order transition between 
the two could be derived this way. 

We argue that this is problematic because the MI tran- 
sition in the actual correlated system is driven by large 
changes in transfer of dynamical spectral weight in re- 



sponse to small changes in the renormalised LDA param- 
eters. In a multi-orbital system, these (the AFE distor- 
tion in VO2) are themselves modified in a priori unknown 
(sometimes non-trivial) ways by electronic correlations. 
Our approach is to search for an instability from the Mott 
insulating (MI) phase (corresponding to one solution of 
the DMFT equations) to the paramagnetic metallic phase 
(the second solution of the DMFT equations) as a func- 
tion of the renormalised AFE distortion, as explained 
above. 

Finally, another source of discrepancy is the use of dif- 
ferent LDA DOS: we have used the ti g DOS as derived 
by Korotin et al. 20 while the LDA+DMFT(QMC) works 
use downfolded LDA DOS. As is clear from a direct com- 
parison between the two LDA results (see Fig. 5), down- 
folding the LDA DOS results in an appreciably larger 
bonding-antibonding feature in the d xy DOS. Clearly, the 
details of the many-body DOS depend sensitively on us- 
ing/not using this device, and the difference between our 
result and those in other works is partly attributable to 
this difference in the corresponding LDA DOS. 
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FIG. 5. LDA partial density of states for the d xy (top) and 
d y z±zx (bottom) orbitals, obtained from Ref. 20 (black) and 
from Ref. 29 (blue). The latter obtained after downfolding. 

While our LDA+DMFT results provide an excellent 
quantitative description of the metallic phase, they give 
only a very limited description of the c-axis dimerized in- 
sulating state (via structural aspects encoded in the LDA 
DOS). This shortcoming has to do with the intrinsic in- 
ability of DMFT itself to treat dynamical correlations 
associated with the quasi-one dimensional spin-Peierls 
transition seen in the I-phase of VO2. Thus, a quan- 
titative description of the I-phase requires an extension 
of our approach to include the above effects via cluster- 
DMFT. Such a program is under way. 

We can however make a few qualitative remarks on the 
dimerized phase. In our analysis, the / phase is driven by 
strong correlations in concert with A a p. Given this, the 
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existence of long- and short V-V distances along c with 
concomitant c-axis spin dimerization would follow as a 
consequence 12 of the / phase and would not be responsi- 
ble for driving the MIT (incorporation of A a p into LDA 
does not yield the I phase) . An analysis along the lines of 
Rcf. 12 is possible, but requires more work within our for- 
mulation. The resulting c-axis spin dimerization in the / 
phase should in principle explain the magnetic properties 
in the / phase. This would also have the generic effect 
of increasing the gap (Eq) to values closer to the exper- 
imental value of 0.6 eV. More work is needed to make 
quantitative comparison, and this is left for the future. 



V. CONCLUSION 

In conclusion, in the light of the theoretical constraints 
imposed by various experimental results, we have stud- 
ied the T-driven MIT accompanied by abrupt changes 
in orbital occupations using the ab-initio LDA+DMFT 
scheme. In complete accord with Mott's ideas, 10 the 
MIT is shown to be of the Mott-Hubbard type, driven 
by an increase in the itinerant carrier density. The or- 
bital occupations in the PM phase are found to be in ex- 
cellent semiquantitative agreement with indications from 
recent polarized XAS studies. Furthermore, good quanti- 
tative agreement with thermodynamic data and the one- 
electron spectrum in both the / and M phases is obtained 
in the same picture, providing a unified description of 
the physics of VO2 across the MIT. Such a description 
can have a wider application to other systems 7 ' 6 where 
strong multi-orbital correlations in concert with struc- 
tural details drive a first-order MIT under external per- 
turbations. 
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